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ABSTRACT 

Rotation curves constrain a galaxy's underlying mass density profile, under the assumption that 
the observed rotation produces a centripetal force that exactly balances the inward force of gravity. 
However, most rotation curves are measured using emission lines from gas, which can experience 
additional forces due to pressure. In realistic galaxy disks, the gas pressure declines with radius, 
providing additional radial support to the disk. The measured tangential rotation speed will therefore 
tend to lag the true circular velocity of a test particle. The gas pressure is dominated by turbulence, 
and we evaluate its likely amplitude from recent estimates of the gas velocity dispersion and surface 
density. We show that where the amplitude of the rotation curve is comparable to the characteristic 
velocities of the interstellar turbulence, pressure support may lead to underestimates of the mass 
density of the underlying dark matter halo and the inner slope of its density profile. These effects 
may be significant for galaxies with rotation speeds < 75kms~ 1 , but are unlikely to be significant in 
higher mass galaxies. We find that pressure support can be sustained over long timescales, because 
any reduction in support due to the conversion of gas into stars is compensated for by an inward 
flow of gas. However, we point to many uncertainties in assessing the importance of pressure support 
in real or simulated galaxies. Thus, while pressure support may help to alleviate possible tensions 
between rotation curve observations and ACDM on kiloparsec scales, it should not be viewed as a 
definitive solution at this time. 

Subject headings: cosmology: dark matter — ISM: kinematics and dynamics — galaxies: ISM - 
galaxies: kinematics and dynamics - galaxies: dwarf 



1. INTRODUCTION 

Galaxy rotation curves offer some of the strongest 
evidence for dark matter on kiloparsec scales 
(jSofue k, Rubinl I2001L and references therein). This 
evidence rests on the assumption that rotating disk 
galaxies are in equilibrium, such that the outward 
centripal force of the rotating disk exactly balances 
the inward gravitational force from the mass interior 
to the disk. The rotation speed seen in atomic gas 
at large radii therefore can constrain a galaxy's mass. 
The apparent "excess" mass inferred from the rotation 
curve, compared to the much smaller mass inferred 
from observations of the baryons within the galaxy, is 
assumed to be due to dark matter. The dark matter 
furthermore must be distributed with an approximately 
p ex r~ 2 density profile in the outer galaxy, to reproduce 
the flat rotation curve that is observed. 

In the centers of galaxies, constraints on the dark mat- 
ter require a more accurate accounting of the baryonic 
mass, since baryons make up a l arger fraction of the mass 
in a galaxy's inner regions (e .g., lKalnaisll 983: Kcntl l 19861 : 
Ivan Albada fc Sancisilll986l) . Unfortunately, assigning a 
mass to the baryonic component of galaxy is uncertain, 
due to our limited knowledge of the stellar mass-to-light 
ratio, the effects of dust extinction, and the contributions 
of atomic and molecular gas. 

Alternatively, the difficulties in probing inner halos 
can be reduced by studying only galaxies with intrin- 
sically low baryon fractions, such as dwarf and low sur- 
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face brightness galaxies. These galaxies have no cen- 
tral bulges, low baryonic surface densities, and high gas 
fractions, thus minimizing the impact of uncertainties in 
the stellar mass-to-light ratio. Detailed analyses of such 
galaxies tend to show more slowly rising rotation curves 
than in massive spirals. Part of the slow rise is due to 
the reduced contribution of baryons to their in ner den- 
sity profiles (although, see iSwaters et al.l 120091 ) . How- 
ever, the slowly rising inner rotation curve is also indica- 
tive of shallow inner density profiles for the dark matter 
halo (p oc r°, with a < —0.25; see the thorough review by 
Ide Blokll2009l and references therein) . These inner slopes 
are in apparent conflict with predictions from numerical 
simulations of dark matter halos in ACDM, which fa- 
vor inner power law slopes with a steeper than -0-75 on 
kiloparse c scales and below (iDubinski fc Carlberg||1991 
Navarro et al 



1997; M oore et al.l 1998: Colin et al.ll2004J: 
Navarro et al.l 12004: lHayashi et a l. 2004; iNavarro et all 
20081 : iStadel et al.1 fcooa K When observational uncer- 
tainties in deriving rotation curves from long-slit spectra 
are taken into account , the appare nt conflict is weaker 
dSpekkens et al.l [200 5: Rhc e et al.l 12004 lHavashi et all 
2004; Valcnzucla et al. 2007), although the discrepancies 
between observations and theory still appe ar to persist 
when full velocity fields are conside r ed (e. g .JSpano et al. 



2001 iKuzio de Narav et"aD [200l [20091 : Ide Blok et al 



2008a). 



In addition to the possible discrepancy in the inner 
slope of the density profile, the overall density likewise 
appears to be significantly lower than seen in simu- 
lated halos (jAlam et al.ll2002T ) . The tension between the 
observed and predicted halo densities ha s been some- 
what reduced by the latest WMAP results (jMaccio et al.1 
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l2008f). but signifi cant offsets remain (see the discussion 
in lSellwoodl[2009l) . 

Taken together, the rotation curve analyses are some of 
the strongest lines of evidence that the success of ACDM 
on large scales breaks down when observed on kiloparsec 
scales. The majority of these analyses rely on the as- 
sumption that the galaxy's radial equilibrium is entirely 
governed by the competition between gravity and cen- 
tripetal force. However, most rotation speeds are mea- 
sured using Hi or Hn emission lines. These lines are emit- 
ted by gaseous disks, which can also experience radial 
forces due to pressure gradients. These pressure gradi- 
ents have the potential to provide an additional outward 
force, which helps to support the disk radially. This sup- 
port allows the disk to remain in equilibrium while rotat- 
ing more slowly than the true circular velocity associated 
with the gravitational potential. Thus, if pressure sup- 
port is neglected, then the inferred dynamical mass of 
the galaxy will be less than the true mass. 

Pressure support in the ISM has several possible 
sources. The most obvious is thermal support (i.e. 
Ptherm = nkT) . However, analysis of the temperatures 
and densities of the multiphase ISM suggest that ther- 
mal pressure falls far short of what is n eeded to support 
the o bserved vertical structure (e.g., iBoulares fc Coxl 
1990). The second potential source of support is mag- 
netic. However, it too appears to be insufficient t o 
support the disk (e.g. Jde Avillez fc Breitschwerdtj [2005). 
The final, and most likely dominant, source of pressure 
su pport comes from interstellar turbulence (see review s 
bv lElmegreen fc Scald [200l iMac Low fc Klessenl 120041 ). 
This turbulence is most likely continuously driven by 
stellar winds and supernovae, although the exact source 
is under some debate, and other driving mechanisms may 
take over in regions where the star formation rate is low. 
The turbulent pressure scales with the gas density and 
the velocity dispersion as Pturb ~ po% , such that regions 
of high gas density and high velocity dispersion have 
higher pressures. Since both of these quantities tend to 
decline with radius in galaxy disks, pressure decreases 
outwards, producing a pressure gradient that can con- 
tribute to the radial support of the disk. 

Corrections for pressure support have a long and varied 
history in the literature. In the early days of interfero- 
metric 21cm mapping of dwarf galaxies, it was recog- 
nized that the turbulent motions within dwarf galaxies 
were comparable to th e amplitude of the rotation curve 
(e.g.. iTullv et al.lll978f ). Thus, pressure support due to 
turbulent motion had the potential to produce a signif- 
icant difference between the true circular velocity and 
the observed rotational speed. It became standard to 
"correct" the observed rotation curve for pressure sup- 
porlo, usually assuming either a gauss ian or a con stant 
radial profile for the turbulent velocity (ITullv et al.1 1978| : 
Skillman et all fl987t ILo et al.l fl993 ICote et all 1200 



Meurer et al.l 119961 ) . These corrections fell out of favor 



2 This correction was often erroneously referred to as a correction 
for "axisymmetric drift", a term that applies to a stellar popula- 
tion's difference from the true circular speed, due to its velocity 
dispersion. The observed offset is produced by the different sign of 
orbital epicycles for stars with guiding centers inside and outside 
of the radius of observation. Although the mathematics is sim- 
ilar, the underlying physical situation is quite different from the 
pressure support of the gas. 



during the last decade, presumedly as the focus turned to 
slightly higher mass low surface brightness galaxies. In 
addition, assessment of the importance of pressure using 
more recent, higher quality data suggested that pressure 
supp ort could be n eglected in most circumstances (e.g., 
ISwaters et~aTll2009D . 

Renewed attention on the importance of pressure sup- 
port comes from the latest generation o f hydrodynamical 
simulations of galaxies. Most notably, Vale nzuela et al.l 
(2007) explicitly measured the difference between the ro- 
tational velocity of the gas and the true underlying cir- 
cular velocity, for model late-type disk galaxies with a 
maximum circular velocity of V c ~ 70 km s~ 1 . They found 
that the observed rotational velocity indeed underesti- 
mated the true circular velocity, and argued that pressure 
support was responsible for a significa nt portion of the 
difference. Although it is likely that thelValcnzucla et al.l 
(2007) simulations were not able to fully resolve all the 
relevant subg rid physics involved in gen erating turbulent 
pressure (e.g. JJoung fc Mac Lowll2"006f ) . their results cer- 
tainly motivate a renewed examination of the effect of 
pressure in disk support, par ticularly in li ght of new ob- 
servational constraints (e.g.. rTamburro et al.l [2009h . In 
addition, recent simulations suggest the possibility that 
estimates of pressure support based on measurements of 
the vertical velocity dispersion may underestimate the 
true pressure supp ort by at least a factor of two (e.g., 
lAgertz et ai]|200g) . 

In this paper we reassess the likely importance of pres- 
sure in different regimes of galaxy mass, for a range of 
plausible galaxy models, using modern data to constrain 
the distribution and kinematics of the gas. We first ad- 
dress the issue analytically, by calculating the expected 
differences between the observed rotation curves and true 
circular velocities, for both NFW halo profiles and re- 
alistic slowly-rising rotation curves. We show that the 
dark matter density profiles inferred from the observed 
rotation curves are extremeley sensitive to the details of 
the assumed pressure distribution. Thus, small changes 
in the assumptions can drastically change the scientific 
conclusions. We then consider whether pressure support 
is a sustainable phenomena. Gas consumption will tend 
to preferentially reduce the gas density in the center of 
galaxies, reducing the pressure support over time. We 
calculate the evolution in the pressure support with time, 
and find that inward redistribution of gas in response to 
decreased pressure support tends to counteract the cen- 
tral gas consumption, such that pressure support can be 
significant over the lifetime of a star forming disk. We 
next analyze the pressure support in a galaxy from the 
THINGS survey. In this particular case, we find that 
pressure support is not significant, as we expected from 
the overall rotation speed of the galaxy. However, the 
exercise highlights the difficulty in accurately assessing 
the amplitude of pressure support in real galaxies. Fi- 
nally, we discuss the possible limitations of application 
of simplified analytic calculations of pressure support. 

2. PRESSURE SUPPORT FOR ROTATION 

The circular speed of a test particle in a circular orbit 
is defined as 



Vc(r) 



d$(r) 



dr 



(1) 
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for an axisymmetric system, where <£>(r) is the gravita- 
tional potential. If the system is in equilibrium, then the 
inward radial force due to gravity must be balanced by 
an equal outward force. Traditionally, one assumes that 
the outward force is entirely dominated by centripetal 
force due to rotation with a velocity Vg , such that 



V 9 2 {r) GM{<r) 



(2) 



where M(< r) is the mass contained with a radius r, 
assuming that the potential is approximately spheroidal 
and thus that there is no significant gravitational force 
from mass at larger radii. This equation implicitly as- 
sumes that all material at a given radius is moving with 
the circular speed, and that that the observed tangential 
velocity Vg obeys 



V e (r) = V c {r). 



(3) 



However, in most cases Vg(r) (the "rotation curve") is 
measured using gas as a tracer. Unlike ideal test par- 
ticles, the gas (usually neutral Hi or ionized Hn) can 
experience an additional outward force due to pressure. 
This gas pressure can potentially provide additional ra- 
dial support to the rotating gaseous disk. The resulting 
outward force allows the gas in the disk to rotate more 
slowly than the circular velocity, while still maintaining 
a stable circular orbit. In such cases, interpreting the 
observed angular velocity Vg as the circular velocity V c 
would lead one to underestimate the mass within r. 

To estimate the importance of pressure support in 
modifying Vg, we consider a parcel of gas at radius r 
within a rotating disk. We assume that the gas has den- 
sity p within a volume dV — dr dA (where dA is the 
surface area of the volume normal to the radius), giving 
a mass dM = pdrdA. For a stable circular orbit, the 
radial forces must be in balance, giving 



= F, 



grav 



F r , 



F F 



(4) 



where F grav , F rot , and Fp are the forces due to grav- 
ity, centripetal acceleration, and the pressure gradient, 
respectively. 

If we assume that the gravitational potential is dom- 
inated by an approximately spheroidal mass component 
(as one might expect for late type galaxies with low sur- 
face density disks), then the radial gravitational force 



~V r $ is 



"grav UlVl 



GM(<r) 



The force due to centripetal acceleration is 

r 
and the force due to the pressure gradient is 



F P = 



— — drdA 
dr 



-dM 



1 dP{r) 



p(r) dr 



(5) 
(6) 

(7) 
(8) 



Substituting back into the force equilibrium equation, 
cancelling dM, and rearranging yields 



KV) = V?(r)-o?(r) 



r dP{r) 
P(r) dr 



(9) 



To calculate the radial pressure gradient, we assume 
that the gas pressure is dominated by turbulent motions, 
rather than thermal processes. One can thus approxi- 
mate the local gas pressure as 



P = Pgas O 2 



(10) 



where p gas is the density of the gas and a 2 is the 
one-dimensional velocity dispersion. Both the gas den- 
sity and velocity dispersion tend to decline with ra- 
dius, leading to a radial gradient in the pressure, which 
produces an outward force. If one assumes that the 
gas disk has constan t thickness (e.g.. lKregel et alJl2004t 
lBottemaet"aT1ll986f PI. then p(r) oc E(r) where E(r) is 
the mass surface density of the gas disk. Equation [9] 
then becomes 



Vj{r)=V 2 (r) 1 + 



!(r) 



V 2 (r) 



dlna r (r) dlnE(r) 



cflnr 



dlnr 



where ay is the radial component of the gas velocity dis- 
persion ellipsoid. 

Equation [TTJ can be used to calculate the true circular 
velocity V c from the observable quantities Vg, E, and a, 
assuming that one can estimate ay from the HI veloc- 
ity dispersion. We will discuss this latter assumption in 
Section 15.11 below. Note also that the size of the cor- 
rection to V c does not depend on the absolute value of 
the gas surface density, and instead only depends on its 
fractional rate of change with radius. 

Equation [TT] suggests that the observed tangential ve- 
locity can differ significantly from the true circular ve- 
locity in regions where the turbulent velocities are of the 
same order as the rotational velocity. In addition, it 
shows the that strength of the pressure support depends 
on the fractional rate of change in the turbulent veloc- 
ities and the gas surface density with radius. Galaxies 
with rapidly declining values of a r or E will have more 
significant pressure support, for the same a r . 

In most applications the full expressions for a r (r) and 
E(r) can be used in equation 1111 However, it is instruc- 
tive to simplify the expression by assuming that the gas 
surface density and 1-d turbulent velocity fall off radially 
like exponentials with scale lengths h r and h a , respec- 
tively. If so, then the force balance equation becomes 



V c 2 (r) = V 2 (r)+a 2 (r) 



(12) 



where we have defined a pressure scale length hp as 



3 Note that the galaxies in both IKregel et al.l l|2004l ') and 
Bottcrna ct al. (1986) have higher average circular speeds (Vg > 
100 kms -1 ) than the galaxies that are relevant for this paper. For 
lower mass galaxies, firm constraints on the radial variation in the 
gas scale height appear to be unavailable, beyond the fact that 
the gas distribution seems to much th icker (compared to the scal e 
length) than in more massive galaxies ( Rovchowdhurv et al. 2010). 



Dalcanton & Stilp 



h P = 



h„ h r 



2h r 



(13) 



Empirically, the gas density falls far more steeply with 
radius than the velocity dispersion does ( Section 13. 1.1|) . 
Therefore, the r/h r term will typically dominate, such 
that hp ~ h r . 

Alternatively, we can attempt to es timate the p r essure 
directly from the gas surface density. Uoung et al.l ()2009f ) 
use hydrodynamic grid simulations to calculate the tur- 
bulent pressure in stratified gas layers, after adopting a 
power-law relationship between gas density and star for- 
mation rate. They find that the turbulent pressure scales 
with the star formation rate per unit area £ s / r like 



rturb OC ^-> s f 7 



0.66±0.04 



(14) 



which, for a Schmidt law (S a / r oc S 14 : iKennicutt] 1 1998( 1 
gives a power-law relationship between the turbulent 
pressure and the local gas surface density: 



±turb OC 2j 



0.92±0.05 



(15) 



Note that the exponent differs from Joung et al.'s equa- 
tion 20 for the relationship between turbulent pressure 
and the gas density p, due to variations in the gas scale 
height with gas density and star formation rate. With 
equation ll5l we convert equation |9] into a form where the 
strength of the a 2 term depends on gas density alone. 



V 2 (r)=V 2 (r)-a 2 (r) 



0.92 r d£(r) 



£(r) dr 



(16) 



For an exponential gas distribution, this equation reduces 
to 



V?(r) = V?(r)-o*(r) 



0.92 r 



(17) 



which is comparable to equation [T2l if hp rs h r /0.92. For 
comparison, our first derivation of the correction to V c 
(equatio n [j"3|) suggests tha t hp » h r /lA for h a ~ 5/i r 
(Section [3X1] below). The Uoung et al.l (|2009f l prescrip- 
tion should therefore produce changes in Vg that are 
smaller than those in equation Q2] by ^60%, for the 
same velocity dispersio n and gas sur f ace d ensity profile. 
Note however that the Uoung et al.l (2009) prescription 
in equation [14] has not yet been calibrated in the low 
star formation efficiency regime, within which the expo- 
nent on the relationship between gas surface density and 
star formation rate is typically mu ch steeper than that 
adopted for a Schmidt law (e.g.. lBigiel et al. 2008). 

Unless otherwise stated, however, we will retain the 
more general form of equation [TTJ when calculating the 
effects of pressure support. 

3. CORRECTIONS TO THE MASS DENSITY 

The reduction of the tangential velocity has a direct 
impact on the inferred mass density Pinferredif). For an 
arbitrary rotation curve, the density at some radius r can 
be calculated from the rotation curve and its derivative 



p(r) 



1 



4ttG 



— 1 



dlnV c 
dlnr 



(18) 



assuming the potential is roughly spherical. 

The mean mass density within a radius r, (p(< r)), 
can also be derived from the rotation curve 



(/>(< r)) 



4nG 



V c {r) 



(19) 



This quanity has the benefit of being a more robust mea- 
surement of the density, since it depends only on the nor- 
malization of the rotation curve, and not its derivative. 
If no corrections are made for pressure support, then 
the mean density inferred to be within each radius is 
reduced from the true density, since Vg(r) < V c (r). For 
the specific case of exponential scale lengths for the gas 
and velocity dispersion, 



a 2 (r) 



(Pinferred(< r)) = {p t rue{< ?*)) 1 



v 2 (r) 



( 2 °) 
confirming that the net result will be to infer a lower 

mean density than the true mean density (ptrue{< r )}- 

Pressure support will also alter the apparent slope of 

the density profile. The power law slope of the density 

a (i.e., p oc r a such that a — , , np ) can be related to the 

V ' r d Inr ' 

d In V c \ 



logarithmic slope of the rotation curve 7 (i.e., 7 = g hi ,. 
through 



2( 7 -l) + 



1 d7 
1 + 7 d In r 



(21) 



Pressure support will tend to make the apparent rota- 
tion curve shallower than the true rotation curve, reduc- 
ing the value of 7. We therefore expect that neglecting 
pressure support would lead one to infer shallower inner 
density profiles (i.e., smaller values of a), for realistic 
values of j(r). 

The reduction in density will be most significant where 
the circular velocity is small compared to the gas veloc- 
ity dispersion, and will thus be noticible over a larger 
fraction of a galaxy's radii when the virial velocity of the 
galaxy is small. Thus, pressure support is potentially 
most significant in the regimes where the apparent dis- 
crepancies from simulated dark matter halos are largest 
- namely, the inner regions of low mass galaxies. 

3.1. The Amplitude of Pressure Support 

To calculate the amplitude of the effects of pressure 
support, we now calculate changes to the rotation curve 
and the inferred density profile. We do so in two ways - 
by starting with an assumed true V c and then calculating 
the rotation curve Vg that would be observed, or by start- 
ing with observations of Vg and then inferring the true V c . 
We consider three baseline models, with v max = 30, 50, 
and lOOkms -1 . For each of the models, we first must 
set the physical parameters of the model galaxies (i.e. 
scale lengths, velocity dispersions, etc). In the following 
sections, we discuss the observations used to constrain 
the scale lengths and amplitudes of the gas density and 
velocity dispersion profiles ( Section 13.1. ip . the circular 
velocities of realistic dark matter halos (Section 13. 1 ,2|) , 
and the observed galaxy rotation curves (Section [3. 1.3|) . 
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3.1.1. Setting the Gas Density and Velocity Dispersion 

Profile 

To assign realistic values for the amplitude of the pres- 
sure gradient, we must first specify the radial behavior 
of the gas density and the velocity disperion profile (e.g., 
Equation [TU|) . We can constrain these values using the 
wealth of Hi observations of ne arby galaxies. 

Based on the observations of iSwaters et al.l (|2002l ). we 
assume that the gas density declines as an exponential 
with a scale length of h r . We choose values for h r us- 
ing observations of the Hi s cale length hn i as a func- 
tion of V c = W20/2 given in ISwaters et a l. (2002) for a 
large sample of dwarf irregular galaxies. We find that 
galaxies follow log l0 (fiHi) = 0.5(14/100 km s" 1 ) —0.1 for 
lOkms - 1 < V c < lOOkms -1 , af ter updating the dis- 
tances in [Swate~& ; Balcells) (|2002[ ) with more recent de- 
terminations derived from the tip of the red giant branch, 
and eliminating galaxies with D > 10 Mpc, for which dis- 
tances are more uncertain. The scatter around the mean 
relation is ah„, =0.12. We note that the Hi scale lengths 
measured by ISwaters et al.l (2002) typically ignored the 
very inner regions, where the Hi profile rolls over to a 
nearly constant surface density of 6 — 10 M© pc~ 2 fol- 
iate- type disks. However, we consider it likely that the 
saturation of the Hi surface density i s due to the con- 
version of Hi into H2 (see Figure 8 of IBigiel et al.l [20081 
for example), in which case the gas density should in 
fact continue rising to the center. We thus assume that 
h r ~ hjil- 

Our next step is to set t he radial veloc ity dispersion 
profile of the galaxies using iTamburro et alJ (l2009f)'s re- 
cent analysis of Hi data from THINGS (jWalter et aTl 
I2008T ). We take three galaxies (from their sample of 11) 
as being representative of the different vel ocity regimes 
we ar e probing - Holmberg II (V max ~ 35; iPuche et al.l 
[1992), NGC 4214 {V max ~ 60 km s~ 1 ; |AUsopdll979l) , an d 
NGC 7793 (V max ~ QOkms-^ ICarignan fc Pu"chel[r990l) : 
note that the exact rotation speeds o f these galaxies are 
uncertain, as ITamburro et al.l ((2009) chose nearly face- 
on galaxies for their velocity dispersion analysis. The 
velocity dispersion profiles are well fit with exponentials 
&z( r ) — Czfiexpi—f/ho-), with a z measured approxi- 
mately perpendicular to the plane of the galaxy. For 
our three representative galaxies, the velocity dispersion 
profile of the gas is well fit with a Zy o = 17.2, 18.2, and 
18.2 km s" 1 and K = 7.3,6.0, and 11.5 kpc, for HoII, 
NGC 4214, and NGC 7793, respectively. There is re- 
markably little variation in the central velocity dispersion 
inferred for the range of galaxy masses we considered, 
and only a modest variation of h a . The scale lengths 
of the velocity dispersion profiles are much longer than 
those inferred for the gas surface density. Using the re- 
lationship between rotation speed and HI scale length 
we derived above, we estimate h a /h r ^4 — 6. Go- 
ing forward, we therefore adopt t7 r o = 18kms _1 for all 
model galaxies, and h a /h r ~ 5. Note that we are forced 
to assume that the velocity dispersion is isotropic (i.e. 
&z = o> ) > given that we lack any reliable way to measure 
a r for the gas. 

3.1.2. Pressure Support for a True NFW Circular Velocity 

With the above profiles for the gas density and velocity 
dispersion, we have sufficient information to calculate the 



radial pressure gradient for our model galaxies. We now 
must make assumptions for the rotational velocities of 
the model galaxies. In this section, we assume that there 
is some true underlying circular velocity profile known 
from N-body simulations, and then calculate the appar- 
ent rotation curve, including the effects of pressure. 

We adopt a true V c (r) for the model galaxies by as- 
suming that the mass distribution is dominated by an 
underlying NFW density profile (Nav arro et al.l Il996f) 
with a virial concentration parameter c V i r . We adopt 
virial velocities V V i r that produce the desired fiducial 
rotation speeds at ~ 10h r (such that typically V max ~ 
1. 25V^ r ). We then use the numerical simulation results 
of IMaccio et al.1 (|2008| ) to fix c V i r as a function of V v i r , 
adopting log 10 (<w) = 1.579 — 0.283 log 10 (14i r /kms~ 1 ) 
for the WMAP5 cosmology; for the three fiducial rota- 
tion speeds of 30, 50, & 100 kms -1 , we derive c V i r = 15.7, 
13.5, & 11.0, respectively. The choice of V v i r and c V { r 
then fixes the core radius r c for the density profile, where 
f c = r V i r /c V i r and (14^/75 kms -1 ) = (r V i r /151 kpc) for 
an assumed WMAP5 cosmology. We then calculate the 
true circular velocity V c (r) expected for the NFW pro- 
file, and derive the apparent angular velocity Vg(r) using 
equation ITTI to include the effects of pressure support. 

We plot the resulting observed rotation curves in Fig- 
ure [TJ for the three fiducial rotation speeds. At high 
rotation speeds V c 3> a r , the angular momentum sup- 
port is far more important than pressure support, and 
the observed rotation curve is an excellent tracer of the 
true circular velocity. At low rotation speeds, however, 
pressure support becomes significant, such that the gas 
disk can rotate more slowly while still being supported 
against collapse. In such cases, the observed rotation 
speed is significantly less than the true circular velocity. 

We calculate the density profiles that would be inferred 
from the observed rotation curves. We plot both the 
true and the inferred density profiles in Figure [U for the 
same rotation curves shown in Figure [T] As expected, 
the inferred density profile is a good match to the true 
density profile for the models where pressure contributes 
little to the overall support of the disk. At lower ro- 
tation speeds, however, pressure support can lead to a 
significant reduction in the inferred density, by a factor 
of 1.46 in the V max — 50 kms -1 model, and 2.18 in the 
V max = 30 kms -1 model. If we consider the la sc atter in 
the halo concentration from IMaccio et al.l (|2008f ) , these 
reductions increase to factors of 1.96 and 6.98, respec- 
tively, at lower concentrations, producing more slowly 
rising circular velocities, and thus larger radial ranges 
over which the ratio of a/V c is large and pressure sup- 
port is significant. At + la higher concentrations, the 
inferred density is reduced by smaller factors of 1.25 and 
1.50 in the V max = 50 kms -1 and 30 kms -1 models, re- 
spectively. The change in the inferred density is far less 
sensitive to our choice of h a /h r . 

Figure [2] suggests that pressure support can poten- 
tially lead to inferred densities that are significantly lower 
than the true underlying density, particularly for galax- 
ies where the velocity dispersion is a significant fraction 
of the rotation speed. However, for the particular form 
of our model galaxies, there is no significant change in 
the inferred inner slope of the density profile. Both the 
true and the inferred density profiles show a cusped inner 
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density profile with a power-law slope of approximately 
-1. Thus, pressure support does not necessarily create 
the appearance of a cored constant density mass profile, 
at least for our particular choice of velocity dispersion 
profile. However, while the literature suggests a possible 
discrepancy between observations of the inner density 
profile slope ( dlnp/dlnr) and N-body simulations (e.g., 
de Blok 2009), we note that such measurements depend 
on the second- derivative of the observed rotation curve 
(e.g., equation IT8)) . and ar e thus highly dependent on 
non-circular mot ions (e.g., Spckkens & Scllwood 2007; 
Wada et al. 2 0021 ) , resolution, and changes in inclination. 

3.1.3. Pressure Support for an Observed Isothermal Circular 

Velocity 

In the above example, we considered how pressure sup- 
port could alter the rotation curve and density profile 
expected for an NFW profile. In this example, we con- 
sider the true circular velocity that would be inferred for 
a frequently used parameterization of the observed rota- 
tion curve shape, after accounting for turbulent pressure 
support. In particular, we consider the arctangent de- 
scription of the observed rotation curve 



Vg(r) =V t ( - ) arctan(r/r t ), 



(22) 



following the notation in iCourteaul {1997). This param- 
eterization produces a rotation curve with an approxi- 
mately linear rise out to a radius r t , beyond which it rolls 
over to a constant "flat" rotation curve with amplitude 
Vt- The associated density profile has a constant density 
core with p w (3/tt 3 G) (V t /r t ) 2 , which then declines as 



2 beyond r t . 



p oc r 

We fixed va lues of r f as a function of V$ jmax us- 
ing data from ICourteaul (|1997h (including the data in 
Tables 6 and 7 for b oth the ICourteaul (|1997[ ) and 
iMathewson et alj (119921) sam p les of late type spirals), 
and from the [S pckkens et al. (2005) sample of dwarfs. 
The lSpekkens et ah (2005) rotation curves have been fit 
with a slightly different functional form, but one that 
agrees extremely well with equation [22] if their distance 
scal e parameters ar e divid ed by 1.2 (i.e. r t « r ex /1.2 in 
the iSpekkens et al.l (|2005l ) notation) and their velocity 
scale parameters are divided by 0.87 (V t ~ V es /0.87). We 
approximate Vg^ max using half of the Hi full-line velocity 
width W20 for all samples. We use only those galaxies 
with asymptotic ally flat rotation curves (J) = for the 
ICourteaul (|1997l ) fits and /3 < 0.03 for the ISpekkens et al.1 
(2005) fits, where ft is a free parameter that accomodates 
rotation curves that do not become asymptotically flat at 
lar ge radii); restric t ing th e sa mple to this subs et allows 
the lSpekkens et al.l (|2005| ) and ICourteaul (|1997| ) samples 
to be easily merged, in spite of their different parameter- 
izations of the rotation curve shape. 

Qualitatively, r t increases for lower rotation speeds 
over the interval from 120 < Vg :max < 300kms _1 , as 
the effect of the central bulge on the density profile 
becomes systematically less important at lower galaxy 
masses. Below 120kms _1 , however, all galaxies are es- 
sentially bulgeless, and the trend reverses, such that r t 
becomes slightly smaller with decreasing rotation speed. 
Because the bulge is essentially gone, the scaling of 
decreasing halo+disk scale length with decreasing halo 



mass dominates at low masses. Since we are only inter- 
ested in low mass galaxies with rotation speeds less than 
lOOkrns -1 , we fit only to this latter regime, and adopt 
log 10 (r f / kpc) = (0.0 ± 0.25) + 0.3 x {Vg, max /100 kms^ 1 ) 
as a good representation of the range of the data. We 
note, however, that the scatter in this regime is quite 
large. Thus, at fixed maximum baryonic rotation speed 
(i.e., V max « W20/2), there is a significant range of den- 
sity concentrations, such that in some galaxies the rota- 
tion curve reaches its peak value quickly, while in others 
the rotation curve is still rising at its last measured point. 

In Figure [3] we show the fiducial observed arctan ro- 
tation curves (Vg(r)) adopted for the three mass ranges, 
along with the adopted velocity dispersion profiles (dis- 
cussed in Section [3.1.11 above). We also show the true 
rotation curves (V c (r)) that would produce the observed 
rotation curve, given an outward pressure force due to 
gas pressure (e.g., equation II lj). As we saw in Figure [1] 
for the NFW density profile circular velocities, the differ- 
ence between the true circular velocity and the observed 
rotation speed are proportionally larger in lower mass 
galaxies. In addition, the differences between V$(r) and 
V c (r) are similar for a fixed maximum rotation speed, 
whether the models began with a true NFW circular ve- 
locity V c (r) (Figured]), or with an observed arctangent 
rotation curve Vg(r) (Figure [3]). 

Significant differences between the two approaches are 
seen, however, when we examine the inferred density 
profiles (equation 1 18 j) for the artangent rotation curves, 
shown in Figure S] The density profile derived from the 
observed rotation curve would naively be interpreted as a 
modified isothermal density profile with a constant den- 
sity core. However, once the effects of pressure support 
are included, the true density profile is seen to have a 
strong central cusp. The differences between the true 
and the observed density profiles are dramatic, and as 
such, failure to take pressure support into account would 
lead to fundamentally different scientific conclusions. 

The relationship between the true and observed density 
profiles are quite different in the two cases we considered. 
When we started with a highly cusped density profile 
and then calculated the observed rotation curve in the 
presence of pressure support, we found little difference 
between the true NFW density profile and that which 
would have been inferred in the absence of corrections 
for pressure support. Both showed steep inner power 
law cusps with comparable slopes. In contrast, when we 
started with an observed arctangent rotation curve, the 
true density profile that would be derived after correcting 
for pressure support is far more cusped than one would 
infer from the observed rotation curve using traditional 
analyses. 

Taken together, the two classes of models suggest that 
the typical effect of correcting for pressure support will 
be to make the inner density profile steeper, although by 
a small amount for apparent rotation curves that already 
suggest a steep inner cusp. More notably, however, the 
inferred density profile can be incredibly sensitive to the 
details of the pressure support. The differences between 
the observed Vg rotation curves in Figures Q] & [3] arc 
rather subtle, but they yield very different amplitudes for 
the correction needed to derive the true density profile. 
The corrections depend on the ratio of Vg to o~ r , and thus 
we expect an equal level of sensitivity to subtle changes 
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in the velocity dispersion profile. 

This sensitivity is likewise reflected in the varying im- 
portance which past authors have placed on pressure 
support. Slight changes in assumptions for the veloc- 
ity dispersion profile (e.g., flat, gaussian, exponential) 
led to a wide range of conclusions about the importance 
of correcting for pressure. Even more worrisome, this 
sensitivity indicates the difficulty in ever making an ac- 
curate correction for pressure support. Minor deviations 
in Vg due to non-circular motions and varying inclination 
will easily compromise reconstruction of the true density 
profile, although numerical experiments seem to indicate 
that such affects are more likely to lead to overestimates 
of the mass, rather than the u nderestimates pro duced 
by neglecting pressure support (Wa da et al.ll2002t ). On 
the other hand, any correction is likely to push the in- 
ferred density profile closer to the true density profile. 
We discuss further complications due to observational 
limitations below in Section [5. f I 

3.2. Reductions in the Mean Density 

For a spherical mass density, the mean interior density 
within some radius r is proportional to (V c (r)/r) 2 . Thus, 
the mean density inferred from a rotation curve will be 
artificially reduced if the observed rotational velocity is 
less than the true circular velocity, as it is in the presence 
of pressure support. This effect can be seen for both sets 
of models in Figures [5] & HI which show systematic re- 
ductions in the density that would be inferred if pressure 
support were not taken into account. 

The mean interior density of a galaxy has been recently 
characterized by a dimensionless density Ay/ 2 > defined 
at the radius ry/2 where the rotation curve has risen 
to half its maximum value V max (i-e., where V@(rv/2) = 
V max /2). The dimensionless density is defined at the 
mean interior density at this radius divided by the cu r- 
rent critical density of the Universe (Al am et al.1 12002'): 
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This quantity can be easily calculated for observed ro- 
tation curves and for simulated dark matter halos. Ob- 
servationally, it is far more robust to uncertainties than 
calculations of the inner slope of the density profile, since 
it depends only on a single point on the rotation curve, at 
radii where the rotational speed is typically much larger 
than non-circular motions. 

Initial comparisons between the value of Ay/ 2 for sim- 
ulated halos and for low surface brightness galaxies sug- 
gested that the mean interior density of the simulated 
CD M halos was much higher than observed in rea l galax- 
ies (jAlam et al.l l2002t iZentner fc Bullock! I2OO30 . even 
when no correction for the baryonic mass of the galax - 
ies was made. Recent comparisons (Ma ccio et al.l 120081) 
show less tension between the predictions and observa- 
tions, largely due to the reduction in the predicted halo 
density with the currently favored WMAP5 parameters. 
However, the two b arred galaxies model ed in de t ail by 
iWeiner et all (|200lD and iZanmar Sanchez et al.l (|2008l ) 



find values of Ay/ 2 that are still an order of magnitude 
below the WMAP5 CDM prediction (although some of 
the low halo density may be due to angular momen- 
tum e xchange from the bar itself; see review by iSellwoodl 
(|2009l) l. 

In Figure [5] we show the predicted measurements of 
the dimensionless density A v / 2 for our fiducial galaxies, 
showing both the true value of the dimensionless density 
(solid symbols), and the artificially low value that would 
be measured in the presence of pressure support (open 
symbols). For both the true NFW models (circles) and 
the observed arctan rotation curve models (triangles), 
pressure support significantly reduces the apparent di- 
mensionless density for low mass galaxies. It has little 
effect for 100 kms -1 halos, however. Thus, correcting for 
pressure support could reduce any observed discrepancy 
for low mass halos, but would have no signficant impact 
at higher masses. 

3.3. Can Pressure Support be Maintained? 

The pressure support discussed above depends on the 
surface density of the gaseous disk. However, star forma- 
tion will tend to erode the gas surface density, and thus 
reduce the pressure support. The gas will preferentially 
be depleted in the inner regions, due to the non-linear 
dependence of the star formation rate on gas density. 
The consumption of gas will then lead to an inward flow 
of gas, which in turn will help restore the pressure sup- 
port. In this section we evaluate the importance of these 
effects, to assess whether pressure support can be a long- 
lived phenomena in an evolving star forming disk. 

The evolution of pressure support in a disk depends 
to first order on the evolution of the gas surface den- 
sity Tig. The gas density at a radius R in a rotating 
disk is governed by a series of differential equations gov- 
erning the flow of matter through the annuli. The first 
of these equations describes the conservation of matter 
during flows: 



dT g 1 d{RV r T g ) _ gS* 
dt + R dR ~ ~ dt '■ 



(25) 



where the gas surface density changes in response to the 
flow of matter between annuli (the second term on the 
LHS) and the conversion of gas into stars, as measured 
by the change in the stellar surface density £*. 

The second equation describes the evolution of the ra- 
dial forces, 
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which allows small radial velocities to develop when the 
angular velocity Vg and pressure P are no longer suffi- 
cient to balance the force that results from the gravita- 
tional potential $. We have neglected the radial advec- 
tion of V r , as terms proportional to V 2 are small. 

The third equation describes constraints from angular 
momentum conservation 



d(RVgTg) 1_ d{RV r ■ RVgT g) 
R 



RVg 



as, 



(27) 



dt R dR dt 

where the change in the angular momentum of the gas 
(first term) is due to angular momentum flux through 
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the annulus at a rate controlled by the radial velocity V r 
(second term) , and to the rate at which angular momen- 
tum drops out of the gas phase and into the stellar phase 
due to star formation (right hand side) . I n this equation 
we hav e ignored the effects of viscosity; see lFirmani et all 
(1996) for a full set of evolutionary equations containing 
viscous terms. Note that a gas disk will tend to collapse 
inwards in the presence of star formation (eqn. I2f3[) , due 
to decreased pressure support resulting from the conver- 
sion of gas into stars (eqn. I25J) . However, the resulting 
inward flow of high angular momentum gas (eqn. I27p will 
increase the circular speed, and thus will help counteract 
the shrinking of the disk. 

To solve this system of equations, we make a series of 
simplifications. We first assume that the star formation 
rate surface density follows a Schmidt-type law, and is 
thus proportional to a power asF of the gas surface den- 
sity, such that <9S*/<9£ = 7/> s / r (S s /So) asF where So is 
a fiducial surface mass density and 'ipsfr is the fiducial 
star formation rate surface density when S g = So- We 
adopt S = 10M Q /pc 2 , ip s f r = 4.07M Q /Gyr/pc 2 , and 
asF = 1-4, in agre ement with the mean trends found by 
IBigielet all (pOOl . 

For our second assumptions, we assume that the grav- 
itational potential is unchanged by the evolution in S 9 
and S*, since the potential will be dominated by the 
dark matter and the stellar disk. We also assume that 
the gas disk starts in gravitational equilibrium (such that 
dV r ]dt = initially) before gas consumption drops the 
pressure support and takes the gas out of equilibrium. 

Third, we assume that because h r /h a ~ 5, radial 
variations in the pressure are driven primarily by the 
change in gas density with radius. We thus assume that 
the velocity dispersion varies slowly enough with radius 
that it can be approximated as a constant, such that 
(l/p g )dP/dR*(a 2 r /Z g )dX g /dR. 

Finally, we assume that the velocity dispersion is con- 
stant with time, and thus that the energy in the turbu- 
lence is continually replenished due to energy input from 
evolving s tars and magneto-rotational instabilities (see 
revie ws bv lMac Low fc Klcssen 2004; Elmegreen fc Scald 
2004). To first order, this assumption has empirical sup- 
port from observations that the average velocity disper- 
sion depends very weakly on star formation rate, varying 
by less than a factor of two when the average star for- 
mation rate per unit ar ea varies by mor e than a factor of 
1000 (c.f., Figure 2 of IDib et al.|[200l . for star forma- 
tion rates comparable to those seen in late type galaxies. 
However, the assumption of constant velocity dispersion 
may be less valid in regions with star formation rates 
higher than > 10 -3 M Q yr _1 kpc~ 2 , for which the appar- 
ent velocity dispersion does appe ar to correlate with sta r 
formation rate (c.f., Figure 1 of iTamburro et al.l f2009). 
although again with only a factor of two range in velocity 
dispersion. Thus, although we expect that the decline in 
star formation rate due to gas consumption could poten- 
tially reduce the turbulent velocity dispersion, and thus 
the degree of pressure support, we expect the temporal 
variation in a to be modest, due to the very weak empiri- 
cal correlations between velocity dispersion and all other 
physical quantities. 

With the above assumptions, we can rewrite the sys- 
tem of equations as: 
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where a ff ,o is the initial gravitational acceleration, set to 
keep the system in balance at the initial timestep (i.e. 
such that dV r /dt = at t — Q). These form a system of 
three differential equations with three unknowns (S 5 (£), 
Vg(t), and V r (tj), which can be integrated forward from 
some assumed initial condition. We adopt an exponential 
profile for S s (t — 0) , with scale length h r chosen accord- 
ing to V m axi as described for our fiducial disks in M3. 1 . It 
We parameterize the initial rotation curve with an arct- 
angent (eq. [22]) , with a scale length chosen as described 
in H3.1.3I We assume that the disk is initially in rota- 
tional equilibrium (V r (£ = 0) = 0), and fix the radial gas 
velocity dispersion at ay — 15 kms -1 . Although it is un- 
likely that a real galaxy would ever be found in such an 
idealized state, these initial conditions are adequate for 
determining the timescales over which pressure support 
is significant. 

In Figure [S] we show the resulting evolution of the 
disk's gas surface density (left column) and rotation 
curve (right column), for fiducial rotation speeds of 
Vmax — 30, 50, & lOOkms -1 (top to bottom), at a se- 
ries of timesteps separated by 1 Gyr (light to dark, with 
the darkest line being the final state). In all disks, the 
gas surface density declines with time, with proportion- 
ally larger reductions in the center due to the higher star 
formation efficiency at high gas densities. 

The models also show noticable radial redistribution of 
gas in the lowest mass galaxies, where pressure support is 
significant. In these galaxies, conversion of gas into stars 
leads to a reduction of pressure support. The reduced 
support then leads to an inward flow of gas, causing the 
gas disk to evolve with a nearly self-similar surface den- 
sity distribution. This inward flow of material raises the 
mean angular momentum at each radius, causing the ro- 
tation speed to increase with time. In contrast, galaxies 
with higher rotation speeds have lower degrees of pres- 
sure support, leading to little radial redistribution of gas, 
and no obvious changes in the rotation curve. 

Although the decline in the gas surface density does 
lead to some radial shrinking of the outskirts of pressure- 
supported gas disks, the rate of shrinking is not dramatic. 
The radius at which the outer disk passes through a fixed 
surface density (say, S 5 ~ 0.1 M Q / pc 2 ) falls by only 10% 
over 5 Gyr for the most pressure-supported galaxy, com- 
pared to 6% for the least pressure supported disk. Thus, 
while the reduction in pressure support in slowly rotat- 
ing galaxies increases the rate of shrinking by a factor of 
two compared to rapidly rotating galaxies, the absolute 
amount of shrinking is small in either case. In constrast, 
the shrinking of disks due t o falling pressure support was 
identified by iStinson et al.l (2009) as a significant effect 
in the early evolution of dwarf galaxies. The rates cal- 
culated here are much smaller than found in the simu- 
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lations, however, presumedly due to some combination 
of more rapid gas consumption and a larger degree of 
pressure support in the simulations. 

The shrinking of the disk helps to refuel gas to the 
central regions. However, like the shrinking in the out- 
skirts, the effect is small. The final surface density in the 
V m ax = lOOkms -1 disk is 17% of the initial value, but is 
somewhat higher (21%) in the most pressure supported 
disk (V max = 30kms _1 ) due to the gas that has flowed 
inwards as the pressure support has decreased. 

The overall conclusion from Figure [6] is that pressure 
support can have an impact on the evolution of gaseous 
disks, but the overall size of the impact should be modest. 
Moreover, while the degree of pressure support does de- 
cline somewhat with time, the overall decline is small, as 
can be seen from the fact that the ratio of cr r /V$ remains 
close to its initial value. Thus, a disk that is partially 
supported by turbulent pressure is likely to be so over 
its lifetime, barring any significant change in its gas dis- 
tribution or star formation rate in response to external 
forces (tidal interactions, infall, etc). 

One caveat to the above conclusion is that the tur- 
bulent velocity of the gas may be time-dependent, as 
opposed to constant as we assumed above. For exam- 
ple, if the turbulent velocity of the gas is proportiona l 
to the star formation rate (e.g., iTamburro et al.l [2009). 
then pressure support should increase during times of 
high star formation activity, and decrease during quies- 
cent intervals. Such a variation could produce epochs of 
disk expansion and contraction, correlated with periods 
of high and low star formation rates, respectively. Empir- 
ically, however, the turbulent velocity does not seem to 
vary by more than a factor of two from galaxy-to-galaxy, 
even when the star formation rate per unit area varies by 
factors of ten or more. Self- regulation may therefore limit 
the strength of any possible correlation between turbu- 
lent velocity and star formation rate, which would then 
keep the degree of pressure support roughly constant. 

4. APPLICATION TO REAL GALAXIES 

The calculations above adopt idealized model galax- 
ies. We now derive the amount of pressure versus rota- 
tional support for real galaxies, using publicly-available 
Hi data. 

Estimating the amplitude of pressure support requires 
measurements of the rotation speed, turbulent velocity, 
and gas surface density. Unfortunately, these quantities 
are frequently difficult to constrain observationally, es- 
pecially within a single galaxy. First, the rotation speed 
is best measured for highly inclined galaxies with well- 
determined inclinations, but the velocity dispersion and 
gas surface density are most reliably measured for face-on 
systems. Second, the radial velocity dispersion controls 
radial pressure support, but only the vertical velocity 
dispersion can be easily measured (although the radial 
component can potentially be derived from detailed mod- 
elling) . Third, the balance between pressure and rotation 
support can be most easily evaluated for velocity fields 
with negligible non-circular motions and constant incli- 
nation with radius, but such systems are rare among low 
mass galaxies. Third, the measured gas surface density 
should include both the molecular and atomic compo- 
nents, but molecular gas is rarely detected in low mass 
galaxies, possibly due to a strong mass-dependent varia- 



tion in the ratio of CO to H 2 (e.g., iBoselli et all 120021) . 

We select galaxies for testing the amplitude of pres- 
sure support from the THINGS sample of high- resolution 
VL A observations of H i (The Hi Nearby Galaxy Sur- 
vey; |WteEiraD[2008). The data cubes for THINGS 
galaxies are publicly available, and rotation curves 
(dc Blok et al. 2008a|) and velocity dispersion profiles 
( Tamburro et al.ll2009f ) are available for the majority of 
the galaxies. We restrict our initial selection to galaxies 
with W20 < lOOkms -1 , since only galaxies with V max ^ 
50kms _1 are likely to show significant pressure support. 
The candidates are M81dwA, Hoi, DD053, M81dwB, 
HoII, NGC 628, NGC 4214, DD0154, and NGC 2366, in 
order of increasing apparent rotation speed (uncorrected 
for inclination). We further restrict the sample to those 
with inclinations > 40°, for which rotation curves can 
be reliably measured; on ly t hese galaxies have had rota- 
tion curves published in Ide Blok et al.1 (|2008a| ) . The re- 
maining galaxies are DDO 154 (V max ~ SOkms^ 1 ) and 
NGC 2366 (V max ~ eOkms" 1 ), both of which are in 
a regime where the effects of pressure support may be 
marginally detectable. However, of these, NGC 2366 has 
a strongly declining rotation curve, suggestive of a kine- 
matically disturbed galaxy, leaving DDO 154 as the only 
candidate for providing a legitimate analog of the models 
discussed in this paper. The few galaxies with potentially 
lower rotation speeds are nearly face-on, have asymmet- 
ric velocity fields, or are nearly spatially unresolved. 

The rotation curve of DDO 154 was kindly given 
to us in electron i c form at by Erwin de Blok, based 
on Ide Blok et al.l (l2008al). W e maintained consistency 
with the Ide Blok et al.l (|2008al) rotation curve by adopt- 
ing their values for dynamical centers, inclinations, and 
distances for the sample galaxies, using the natural- 
weighted data sets. 

We estimate the average Hi surface density T,hi and 
velocity dispersion ohi as functions of galactic radius. 
We obtained the inclination-corrected surface density 
from the moment map, and the velocity dispersion from 
the moment 2 map, by integra ting elliptical annuli w ith 
a width equal to that used by Ide Blok et al.l (|2008a| ) to 
derive the rotation curve. The total gas surface density 
at each radius was assumed to be T, gas (r) = 1.4xEjj/(r), 
where the factor of 1.4 is an approximate correction for 
metals and molecular gas (since the latter is currently 
unavailable for these galaxies). We also considered a 
more extreme correction for molecular gas by assuming 
that T. gas follows an exponential profile into the center 
of the galaxy, such that the outer regions are dominated 
by atomic gas, but the inner regions are dominated by 
molecular gas, producing the apparent roll-over of Y>hi 
towards the center. Both of the resulting gas density 
profiles are shown in Figure [3 

The pressure at each radius was taken to be P(r) = 
(cr| f/ (r)E ffas (r)), which implicitly assumes that a r s=s 
a hi ■ This assumption is not necessarily true, given that 
interaction between the ISM and stellar feedback may 
well differ between the vertical direction (where the ISM 
becomes diffuse over 100 pc scales) and the radial and 
tangential directions (where the ISM remains dense on 
kiloparsec scales). This calculation also assumes that the 
pressure is due to random motions of the gas, rather than 
locally coherent motions due to expanding superbubbles. 
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To better capture the broad radial trend in the pressure 
gradient, we chose to smooth the velocity dispersion and 
surface density before calculating the pressure. We then 
use equation[TT]to calculate the true circular velocity that 
would be observed in the absence of pressure support. 

In Figure [7] we plot the resulting data for DDO 154. 
The right hand panel shows the velocity dispersion pro- 
file (dotted line) , the observed tangential velocity (black 
solid line), and the inferred underlying circular velocity 
(red), both for the observed gas surface density (solid 
line) and for an exponential surface density extrapolated 
to the center (dashed lines). There is little evidence 
that pressure support plays any significant role in this 
particular galaxy. The pressure gradient appears to be 
negligible in the inner regions, due to the relatively flat 
velocity dispersion profile and inner gas surface density 
profile. Moreover, the Hi velocity dispersion has a low 
amplitude compared to the fiducial model (~ lOkms -1 
versus 18kms _1 ). There is some evidence for pressure 
support in the outer regions, for which the true rotation 
speed may be ~10% higher than observed. However, 
we have verified that this difference makes no significant 
changes in the inferred density profile. Likewise, the in- 
ferred dimensionless density Ay/2 increases by less than 
2%, to log 10 Ay/ 2 = 5.26, which still lies below the fidu- 
cial curve from simulations in Figure [5] However, while 
the evidence for pressure support is small in this case, 
this galaxy is not in a regime where we expected pres- 
sure support to be obvious. 

5. DISCUSSION OF POSSIBLE LIMITATIONS 

Our analysis of pressure support reveals a number of 
contradictions. On the one hand, straight-forward an- 
alytic calculations demonstrate that pressure support is 
likely to contribute significantly to the radial support 
of low mass galaxy disks, in agreement with numerical 
simulations. On the other hand, observational evidence 
for pressure support is weak. Our analysis of DDO 154 
finds that corrections for pressure support would make 
only minor chan ges to the interpretat ion of the observed 
rotation curve. ISwaters et al.l (|2009l ) rea ches a similar 
conclu sion for the dwarf irregulars in the ISwaters et al.l 
(2002) sample, which have 14>60kms _1 . 

These apparent contradictions must arise from limita- 
tions in the various analyses. Each of the approaches 
for assessing pressure support come with a number of 
caveats, which we now discuss. 

5.1. Limitations of Observations of Pressure Support 

Assessing the importance of turbulent pressure support 
requires knowledge of the radial turbulent velocity dis- 
persion, the gas surface density, and the tangential rota- 
tion speed. Unfortunately, observations cannot measure 
all of these quantities within a single galaxy. Instead, 
we are forced to make several approximations, none of 
which may be completely satisfactory. 

The first approximation is the use of the line-of-sight Hi 
velocity dispersion as a substitute for the radial velocity 
dispersion. There is no reason to believe that the ve- 
locity dispersion of the gas is strictly isotropic, given the 
very different boundary conditions for motions within the 
plane compared to motions in the vertical direction. The 
amplitude of o~ r may therefore be systematically differ- 
ent from the line-of-sight dispersion measured in nearly 



face-on galaxies. Adaptive m esh simulations of rotating 
disks by lAgertz et al.l {2009) find that the velocity dis- 
persion in the plane of the galaxy is at least a factor of 
two higher than the vertical component of the velocity 
dispersion. Because the importance of pressure support 
scales as o~ r /V c , using the measured vertical velocity dis- 
persion (rs a z ) may significantly underestimate the de- 
gree of pressure support, as can be seen from the pref- 
actor of the term in brackets in equation [9] Assuming a 
factor of two difference between a r and o~ z , the degree of 
pressure support for a galaxy with V c ~ 70kms _1 may 
actually be as high as that calculated for a galaxy with 
14 ^35 km s^ 1 . Note that the term within the brackets 
in equation [9] will be unlikely to change significantly if o~ z 
is used as an estimate for o>, since (1) the logarithmic 
gradient in a is likely to be similar between the two di- 
rections, assuming that the shape of the velocity ellipsoid 
varies slowly with radius; and (2) the gradient in a con- 
tributes less to the pressure gradient than the gradient 
in surface density, and thus uncertainties in the velocity 
dispersion gradient are not likely to dominate the error 
budget. 

Also concerning is that the line-of-sight Hi velocity 
dispersion only tracks one particular phase of the gas. 
Thus, Hi velocity dispersions are likely to be biased trac- 
ers of the velocity dispersion in regions of the galaxy that 
have significant molecular or ionized gas fractions. These 
other phases are known t o have different sca le heights 
than the Hi (see review by Ivan de r Hulst 1996), and thus 
must have different characteristic velocity dispersions. 
Unfortunately, data on these other phases have histori- 
cally been far more difficult to obtain than HI. Some hope 
for using only HI observations com es from recent grid - 
based hydrodynamic simulations by Uoung et al.l {2009) , 
who find that different gas phases tend to come into 
turbulent pressure equilibrium (mirroring what has long 
been known for their thermal pressures). Thus, mea- 
suring the atomic turbulent pressure alone may provide 
an indirect estimate of the pressure from molecular and 
ionized phases. 

In addition to the challenge of measuring turbulent ve- 
locities, there are also difficulties in comparing the ob- 
served tangential velocities to analytic calculations, due 
to the presence of non-circular motions. Analytic cal- 
culations assume that the gas moves on perfect circular 
orbits, and thus has no radial streaming motions which 
could otherwise support the radial extent of the disk. 
However, most late-type galaxies show some evidence 
for non-circu l ar motions at the 5-20 km s -1 level (e.g., 
iGentile et all 120051: pTrachternach eFaLl 120081 ). which is 
comparable to the typical amplitudes of turbulent ve- 
locities. Thus, when significant non-circular motions 
are present, the assessment of the importance of pres- 
sure support becomes more complicated than assumed 
in most calculations. Unfortunately, it seems that non- 
circular motions are almost always present in dwarf 
galaxies, particularly at the low masses where pressure 
support is most likely to be significant. In our search for 
a suitable galaxy to compare with our calculations (Sec- 
tion SJ, it proved to be difficult to find low mass galaxies 
for which the velocity fields do not appear to deviate from 
what is expected for a simple rotating disk. Signatures of 
warps and elliptical streaming motions were sufficiently 
common that only one system barely qualified. 
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The difficulty of finding approrpiate systems for evalu- 
ating pressure support is even worse in lower mass galax- 
ies than DDO 154. Although these systems would poten- 
tially have the most dramatic evidence for pressure sup- 
port j_the^_alsoliaveJri£rea^in£]v_ch^ velocity fields 
(see ILo et al] [l993; Begum ct al. 2008), making the mea- 
surement of the tangential velocities as complicated as 
measuring the radial velocity dispersion. Furthermore, 
the gas disks of low mass galaxies frequently appear to be 
far from equilibrium. The p resence of dramatic H i holes 
(e.g.. see reviews by Brinks ct al. 2007; Ivan der Hulstl 
1996) and complicated velocity fields indicates that the 
entire basis of the calculation (i.e., a disk in equilibrium) 
may be suspect. 

There are equally large limitations in measuring the 
gas surface density. The easiest gas component to mea- 
sure is atomic Hydrogen through the 21cm line. How- 
ever, molecular and ionized gas also contribute to the 
turbulent pressure, but are far more difficult to measure. 
Molecular gas in particular is extremely challenging to 
identify in late-type galaxies, for which there are few CO 
detections, in spite of copious evid ence for star formation 
(lisrael et al.lll99 5: Mat thews et aJJl2005t iDas et all 2006: 
iLeroy et al.1120051 ). The lack of CO detections most likely 
arises from factor of ten variations in the Xco conversion 
factor from CO to H2, rather than any lack of molecu- 
lar gas fe.g..lOhta et al.lll993l:lA7imoto et al.lll996Hlsrae]| 
1997; Ma dden et al.ll997HBoselli et al.ll2002l : LLerov et al.l 
2007). Thus, the true contribution of H2 to the turbu- 
lent pressure is particularly difficult to constrain in the 
systems where it is most important. Approximating the 
molecular gas distribution by scaling the Hi observations 
is unsatisfactory as well, since the ratio of H2 to Hi is 
known to vary spatially within galaxies (jWong fc Blitd 
12001 ILeroy etaLl[20Qg , e.g..). 

The ionized gas component is as difficult to constrain 
as the molecular. Hot SN-heated gas is likely to fill the 
Hi "holes" that are prevalent in dwarf galaxies. However, 
x-ray detections of such gas are prohibitively expensive. 
On the other hand, hot gas appears to be a small fraction 
of th e total mass in the Milky Way disk (e.g., iFerrierd 
I1998T ). suggesting that it may not necessarily be a sig- 
nificant source of pressure in dwarfs, provided one can 
safely extrapolate these results to such a different mass 
regime. 

5.2. Limitations of Simulations of Pressure Support 

The pressure support in a stable galaxy disk depends 
on the properties of the interstellar turbulence. This 
turbulence is thought to be driven by various forms of 
stellar feedback (winds, SNe, etc), which inject energy 
int o the ISM over a range of physical scales (see review s 
bv lElmegreen fc Scald l200l iMac Low fc Klesser3l200l . 
This energy then cascades to an even wide range of scales, 
forming a broad power spectrum of turbulent energies. 
The pressure that results from the turbulence thus de- 
pends on the sum of kinetic energy across a wide range 
of length scales. 

Unfortunately for simulations, it appears that the vast 
majority of the kinetic energy is d ue to tubulent mo- 
tions on scales smaller than ~200 pc ([Joung fc Mac Low! 
2006). These small scales have been successfully re- 
solved using either adaptive mesh refinement (AMR) 
or smoothed particle hydrodynamic (SPH) codes (e.g., 



iPrice fc Federrathl 120101 ) . but only for simulations that 
model a small portion of the disk. While these high- 
resolution simulations can resolve the physical scales that 
dominate the turbulent energy, the same codes have diffi- 
culty simultaneoulsy simulating the large scale behavior 
of the galactic disk and its cosmological environment. 
In constrast, simulations that capture the physics of the 
disk on larger scales cannot currently fully resolve the 
"subgrid physics" of turbulence. 

To explore the suitability of simulations for assessing 
the amplitude of pressure support, we will consider two 
state of the art sim ulations. The first is the suite of SPH 
simulations used by I Valenzuela et ail (|2007[ ) to argue for 
the significance of pressure support, feedback, projec- 
tion effects, and non-circular motions in explaining the 
discrepancy between observed rotation curves and the 
canoni cal NFW profile. T he second is the AMR simula- 
tion of Uoung et al.l (|2009| ) used to assess the scaling be- 
tween pressure and star formation rate in self-gravitating 
stratifi ed gas layers. 

The iValenzuela et al.l (|2007l) SPH simulations mod- 
eled a moderately low mass dwarf galaxy, created to 
be an analog of nearby dwarf irregular galaxies (V c ~ 
70kms~ 1 ). Their initial conditions embedded a mixed 
gas+star disk within a live axisymmetric NFW dark mat- 
ter halo. Their m odel was then evolved using the SPH 
code GASOLINE (jWadslev et all l200l for 1 Gyr. Be- 
cause their simulation considered an isolated disk evolved 
for a short period of time, the galaxy is unaffected by 
any larger scale gravitational forces or long term evolu- 
tion driven by interactions with other galaxies. These 
restrictions allowed them to use a large number of parti- 
cles (2 million) and a small softening length (30 pc), such 
that Newtonian forces were reached beyond 60 pc. The 
resulting gas disk successfully reproduces the holes and 
filaments that are characteristic of many dwarf galaxies, 
and has a slowly rising rotation curve, in spite of being 
embedded within a true NFW profile. 

Unfortunately, even these carefully crafted simulations 
do not seem to be fully capturing the effects of pressure 
support. The midplane velocity field (their Figure 1) 
shows coherent rotational flow with significant radial ve- 
locities in the inner regions. However, it shows none of 
the fully turbulent motion that likely dominates the in- 
terstellar pressure. Comparing their effective force reso- 
lution (~ 60 pc) to the scales where the turbulent kinetic 
energy dominates (< 200 pc) suggests that they are un- 
likely to have resolved the full amplitude of the turbulent 
pressure!. Thus, their success in reproducing slowly ris- 
ing rotation curves is likely due to the dominance of the 
other effects they identified, rather than turbulent pres- 
sure. While this work leaves open the question of the 
amplitude of pressure support, it does point to mecha- 
nisms which may well cause the apparent discrepancies 
in higher mass galaxies where pressure is expected to be 
negligible. 

In constrast, the Uoung et al.l (|2009T ) AMR simulations 
do a superb job of resolving the turbulent motions down 
to very small scales (~2pc). However, this high spa- 

4 Note that the difficulty in resolving the pressure is likely to be- 
come even more severe at higher star formation rates, which appear 
to shift the ki netic energy spect rum to even shorter characteristic 
length scales (Joung ct al. 2009). 
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tial resolution comes at the expense of not modeling the 
full galactic disk. Instead, they simulate small patches 
of a vertically stratified gas layer, 0.5 kpc on a side, with 
periodic boundary conditions in the plane. These simula- 
tions therefore are unable to capture radial flows or shear 
fields, and must view the disk in steady state. There are 
also several limitations in the simulation, such as the lack 
of correlation between local density and the local SN rate, 
which may contribute to the simulations' apparent fail- 
ure to reproduce the large holes seen in HI observations. 
On the other hand, the lJoung e t al. (2009) simulations 
do offer a number of lessons for thinking about turbulent 
pressure in general. First, both the thermal and turbu- 
lent pressures appear to be in pressure equilibrium across 
phases, such that the turbulent pressure of the hot gas 
is comparable to the turbulent pressure of the cold gas 
on the larger scales where the turbulent pressure domi- 
nates. The net result is that pressure shows much lower 
variation with position than either density or temper- 
ature. Second, the turbulent pressure scales with the 
star formation rate, since both are correlated with en- 
ergy input from SNe. This scaling relationship provides 
a means to potentially capture the subgrid physics of 
turbule nce, which is other wise unresolved in SPH simu- 
lations ()Joung et al.ll2009h . The next generation of SPH 
simulations may thus have the power to capture both 
the large scale galactic context provided by SPH, while 
better assessing the degree of pressure support. 

5.3. Limitations of Analytic Calculations of Pressure 

Support 

The analytic calculations presented in Section[2]require 
a number of compromises that are made for analytic ex- 
pediency, but that probably fail to capture some of the 
relevant physics. First, as discussed in SectionEj] gas in 
real low mass galaxies shows significant deviations from 
pure circular motion, in contrast to what is assumed in 
Section [2] Bulk non-circular motions add an additional 
term to the Jeans equation, which is not accounted for 
in our calculations. The simulations in I Valenzuela et al.l 
(2007) show a significant contribution from non-circular 
motion, that leads to a reduction in the inner rotation 
curve. Second, the presence of holes and warps suggest 
that low mass galaxy disks are frequently not in a steady- 
state equilibrium, again in conflict with the analytic as- 
sumptions. Third, the derivations in Sections [2] & [3] as- 
sume that the gas disk has no significant viscosity, ve- 
locity anisotropy, or tilt in the velocity ellipsoid. Finally, 
we ignore any support due to thermal pressure, which 
dominates only on smaller scales and lower pressures 

6. CONCLUSIONS 



The calculations above suggest that in low mass galax- 
ies with V c < SOkms" 1 , gas pressure may help to sup- 
port the disk against the inward force of gravity. In such 
cases, the observed tangential velocity of rotation curve 
tracers like Hn or Hi will be smaller than the circular 
velocity of an idealized test particle, leading to under- 
estimates of the mass contained within a given radius. 
Accounting for pressure gradients in the gas may lead 
to higher mean densities and steeper inner halo profile 
slopes for low mass galaxies, thus reducing possible con- 
flicts with ACMD on small scales. Pressure support can 
be maintained over long time scales, even in the pres- 
ence of star formation. In such cases, the reduction in 
pressure support due to central star formation leads to 
an inward flow of gas, which then re-establishes pressure 
support. 

However, while our calculations suggest that pressure 
gradients could help support the centers of low mass 
galaxies, we find that there are many obstacles to ever 
establishing convincing observational or numerical evi- 
dence for pressure's importance. Pressure support de- 
pends on the radial velocity dispersion of the entire ISM, 
but measurements are typically limited to measuring 
only the vertical velocity dispersion of the Hi. In ad- 
dition, the effects of complex non-circular motions and 
evolving superbubbles may dominate the structure of the 
gas disk, making corrections based on simple steady- 
state circular models invalid in the low mass galaxies 
within which we expect pressure to have the largest im- 
pact. Likewise, simulations currently do not have suffi- 
cient resolution to simultaneously capture both the large 
scale galaxy dynamics and the small scale interstellar 
turbulence which dominates the pressure. In light of 
these caveats, while we consider pressure support to be 
an intriguing mechanism to reduce tensions with ACDM 
on small scales, it is by no means a conclusive solution 
to current concerns. Other proposed mechanisms (e.g., 
iGovernato et al.lfeOlOl ) may play an equally or more im- 
portant role. 
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Fig. 1. — True (heavy black; V c ) and apparent (light red; Vg) circular speed, including the effects of pressure support, for a fiducial 
intrinsic NFW density profile in the WMAP5 cosmology, and assuming h a /h r = 5. The velocity dispersion profile is plotted with a dashed 
line. The difference between the true and observed rotation speed is more pronounced for less concentrated halos (which have more slowly 
rising rotation curves), and for larger values of h„-/h r (not shown). 
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Fig. 2. — True (heavy black; ptrue) and inferred (light red; pi„f erred) density profile, including the effects of pressure support, for a 
fiducial intrinsic NFW density profile in a WMAP5 cosmology, and assuming h a /hr = 5. The difference between the true and observed 
density is more pronounced for less concentrated halos, and for larger values of h a jh r (not shown). 
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Fig. 3. — True (heavy black; V c ) and apparent (light red; Vg) circular speed, including the effects of pressure support, for an apparent 
arctan rotation curve, assuming h a /h r = 5. The velocity dispersion profile is plotted with a dashed line. The difference between the true 
and observed rotation speed is more pronounced for larger values of h„ /h r , and for larger values of rt at fixed V m ax (not shown) . 
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Fig. 4. — True (heavy black; ptrue) and inferred (light red; pi n f erred) density profile, including the effects of pressure support, assuming 
hc/hr = 5; for this case, the observed rotation curve was fixed to follow an arctan function, consistent with a constant density core. 
Although the "observed" density profile (i.e. what one would infer from the observed arctan rotation curve) has a flat inner core, the true 
density profile (inferred after corrections for pressure support) has a central cusp, with an inner power law slope of -0.96, -0.7, and -0.35 
for the 30kms _1 , 50kms — 1 , and lOOkms -1 halos, respectively. The inner central slope becomes slightly steeper for larger values of rt or 
of h a /h r . 
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Fig. 5. — Dimensionless density as a function of rotation speed, for observed (open) and true (filled) rotation curves. Circles are for 
rotation curves based on true NFW density profiles (Figure [Tj. and triangles are fo r observed arctangent rotation curves (Figure [3}. The 
dotted line is the approximate relationship for WMAP5, from Maccio ct al. (2008), for reference. Failure to correct for pressure support 
can lead to a factor of ten underestimate in the densities of low mass dwarf galaxies. In contrast, the inferred densities of higher mass 
galaxies (Vc^ lOOkms -1 ) are unaffected by pressure support for our fiducial models. 
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Fig. 6. — Evolution of gas surface density (left) and rotation speed (right) in rotating disks, includirig the effects of pressure support. 
The gas disk starts with an exponential surface density (light line), and evolves due to gas consumption and radial redistribution in 
response to declining pressure support. Successively darker lines show the evolution in 1 Gyrtimcstcps, for fiducial halos with V m ax = 
30, 50, & lOOkrns - 1 (top to bottom) and a constant turbulent velocity of ay = 15kms _1 . Evolution in the gas density due to star 
formation does not drastically alter the degree of pressure support in the disk, due to the resulting radial inflow of gas from the outer disk. 
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Fig. 7. — Gas surface density (left) and velocity profiles (right) for DDO f54, derived from Hi data from the THINGS survey. The left 



hand panel shows the inferred gas surface density profile, assuming that E s 



1.4 X T,ui- The dashed line shows the gas profile if one 



assumes that the total gas density profile is an exponential, with a center dominated by undetected molecular gas. The right hand panel 
shows the observed Hi line-of-sight velocity dispersion profile (dotted line) , and the observed tangential velocity (solid black line) , corrected 
for inclination. The red lines show the effects of correcting for pressure support assuming an exponential gas distribution (dashed line) 
and the observed gas distribution (solid line). Pressure support appears to be negligible in the inner regions of this particular galaxy, as 
expected based on the models in Figure [4] 



